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Abstract 

Pseudoscalar particles, with almost zero mass and very weak coupling to the visible matter, arise in 
many extensions of the standard model of particle physics. Their mixing with photons in the presence 
of an external magnetic field leads to many interesting astrophysical and cosmological consequences. 
This mixing depends on the medium properties, the momentum of the photon and the background 
magnetic field. Here we give a general treatment of pseudoscalar-photon oscillations in a background 
magnetic field, taking the Faraday term into account. We give predictions valid in all regimes, 
under the assumption that the frequency of the wave is much higher than the plasma frequency of 
the medium. At sufficiently high frequencies, the Faraday effect is negligible and we reproduce the 
standard pseudoscalar-photon mixing phenomenon. However at low frequencies, where Faraday effect 
is important, the mixing formulae are considerably modified. We explicitly compute the contribution 
due to the longitudinal mode of the photon and show that it is negligible. 



1 Introduction 

The standard model of particle physics suggests the existence of a very light, weakly coupled particle, 
called the axion. It arises as a pseudo-Goldstone boson of the broken Peccei Quinn (PQ) symmetry in a 
generalization of the standard model [1-7]. Similar particles are also predicted by supergravity [8] and 
superstrings theory [9,10]. Such a pseudoscalar particle has an effective coupling to two photons. As 
a consequence, in an external magnetic field, an axion can oscillate to a photon and vice versa [11-18]. 
The mixing can change both the intensity and the state of polarization of the photons. Since the mixing 
is dependent on the frequency of the wave, it also leads to a change in the spectrum of electromagnetic 
radiation [19-23]. This phenomenon has been used for laboratory and astrophysical searches for such 
particles, leading to stringent limits [24-34]. Furthermore, the astrophysical and cosmological conse- 
quences of this mixing have been studied extensively in the literature [33-44]. The mixing increases the 
transparency of the intcrgalactic and galactic medium to propagation of high energy photons due to the 
very weak coupling of pseudoscalars to visible matter [45-47] . 

Inside a medium, there exists an additional loop induced, axion-photon vertex, providing an extra 
contribution to the usual axion-photon coupling in vacuum. In some kinematical limit this contribution 
has medium, momentum and magnetic field dependence. In vacuum one of the two transverse modes of 
photon gets coupled to an axion. In a medium the photon acquires an extra degree of freedom (i.e., the 



longitudinal mode) and as a consequence axion also couples to this additional degree of freedom. In a 
magnetized medium, the two transverse degrees of freedom are usually coupled due to the Faraday effect. 
Hence in the presence of axions the two transverse modes would also get coupled to the longitudinal 
degree of freedom. 

In this paper we study the axion-photon evolution, taking the Faraday term into account. Our 
objectives in this study are two-folds. First we give a general treatment of pseudoscalar-photon mixing 
in a medium. This allows us to compute the next to leading order temperature dependent corrections. 
Next we give general solutions to the oscillation problem, taking into account the Faraday effect. At very 
large frequencies we expect that the Faraday effect would be negligible and the standard treatment of 
pseudoscala-photon mixing would apply. However at low frequencies the standard pseudoscalar-photon 
mixing result may deviate significantly due to the presence of near degeneracies in the mixing matrix. In 
several laboratory and astrophysical situations this regime may be applicable and hence our treatment 
would be useful in such cases. In the present paper we restrict ourselves to providing a general treatment 
and do not address the issue of applications to laboratory experiments or to astrophysics. 

1.1 Interactions and the Polarization Tensor in a Medium 

The classical Lagrangian of a free electromagnetic field, including the gauge fixing term is given by, 

C = ~F^(x)F^(x) + 1 (d a A a (x)) 2 +T t .A . (1.1) 

In Eq. (|l.l|) . F^(x) = (9 M A I/ (x) — 9 1/ A A '(x)), £ is the gauge fixing parameter, j® xt is the external current 
and A(x)'s are the vector potential. For the sake of simplicity, we use the Feynman gauge and set £ = 1. 

As one takes quantum corrections into account, the quadratic part of the tree level Lagrangian gets 
modified because of quantum corrections coming from terms proportional to the vacuum polarization 
tensor n^^. The resulting Lagrangian, in the momentum space, is 

£ = \ [-*%* + IV(fc)] A»(k)A"(-k) - j^*A"(fe) + Co. (1.2) 

In Eq. (|1.2[) above, Lq corresponds to the gauge fixing term and g^ = (g M „ — ) . In presence of a 
medium of finite density and temperature, the polarization tensor get corrections from the matter and 
temperature dependent parts. The resulting equation of motion for such a system can be written as, 



[-k 2 ~g^+n^{k)]A v {k)=f^, (1.3) 

where we have not retained the pieces coming from the gauge fixing term (i.e. terms proportional to ^pr 1 ) 
since they can be shown to vanish. In a medium composed of electrons, the polarization tensor, H flu (k), 
can be written in terms of scalar form factors and tensors composed out of the four vectors available 
for the system, i.e., the medium four velocity = (1,0,0,0) and k^ = (uj,k). They are expressed as 
follows [48], 

n^(fc) = U T R^ + n L Q Al „, 



f u 'I u 



where <^ n - ~- n ( L4 ) 

in absence of any external field. The vector it M is given by u M = g^vu" . The scalar form factor n^(fc), 
corresponding to the longitudinal degree of freedom is given by, 

U L (k) = -^-U^(k)u^u\ where, u'VTV(fc) = w* ( ^ + 3^f- ] (1.5) 
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Similarly the transverse form factor Tlx is given by the following expression: 

II T (fc) = R^U^(k) and, R^IV(fc) = [l + i*L — ] . (1.6) 



In the expressions above uj p denotes the plasma frequency. In the classical limit, to leading order in 
it is given by: 



-v-fT'O-s) (L7) 

In Eq. (|1.7|) n e is the number density of electrons. 

We next express the form factors in terms of the dielectric constants. Denoting, F^ v = —i(k^A v — 
/Cj,A M ), so that E = iuiA — ikA° and B = ik x A, one can further define longitudinal and transverse parts 
of the electric field as. E\ = k(k.E), E t = E — E\. Using the Kubo formuki (1) l for linear response analysis, 
the induced current, j md (k), can be written as, 



j lnd = iuj 



(1 - ei)Ei + (1 - e t )Et\ (1.8) 
where longitudinal dielectric function ei and transverse dielectric function e t are given as follows, 

(l-^) = 5 (1-9) 

Using the relations given by Eq. (| 1 .6(1 for IIt(A;) and IlL(fc), one can further express the dielectric 
functions in terms of u) and k. 



2 Inclusion of the Faraday Term 
2.1 Exact Faraday Term 

The contribution to the vacuum polarization tensor which is odd in B has been estimated using real time 
finite temperature field theory in Refs. [49,50]. It can be expressed as, 



KM 



Me 



(2tt)4 
Aie 2 e\ pa ,pk 13 



v-(p) 

d 4 P 



ds e*( p " 



ds' e* (p '' s,) 



+i? (2a) 



(2tt)< 



r]-{p) / dse*M / ds' e*( p '< s '> 



Z „ ,Z „ i t&neBs taneSs' , ~ 
p Q " taneBs + p'"" taneSs' : nl ■ „ (p + p') 11 



tane£?(s + s') 



This is exact to all orders in eB. 



(2.1) 



In Eq. (|2.ip ay stands for or the 3. Moreover appearence of a\\ and together in any product would 
mean: if a\\ takes the value 0, would take the value 3 and vice versa. Furthermore for thermal fermions, 
t]f(p) is the distribution function, 



r) F (p) = 9(p-u)f F (p,(i,0) + 0(-p ■ u)f F (-p, -ft, 0) 



(2.2) 



(i) 



i.e. j^(fc) = -IV(fc)yr(fc). 
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Here, G is the step function, which takes the value +1 for positive values of its argument and vanishes 
for negative values of the argument, and fp denotes the Fermi-Dirac distribution function, 

1 



;/3(p-li-Al) _|_ 1 



(2.3) 



Lastly <f>{p, s) in the exponential stands for, 



, 2 tan(eSs) 2 
$0, s) = is[pf ] ^ p±-mT I 



(2.4) 
(2.5) 

However performing the integrals in Eq. (|2.ip and arriving at a compact form is difficult. This integral, 
can be performed in the long wavelength limit. However in this case, the expression does not look gauge 
invariant unless the explicit long wavelength limit is maintained. 



2.2 Equation of Motion with Faraday Contribution 

In presence of external mgnetic field, Eq. (|1.3[) is further modified to [49,50] 

[-k 2 ~ 9 ^ + n„„(k) + n^(fc)] A v {k) = i-*. 



We define, 



In the limit \K\ — > 0, one can express, 



\K\ = 



n p (fc)P^, 



with 



fcC 
A 



/3 



Cti 



K I— i-O 



The limit |A"| — > should be taken in such a way that 

W\ 

The scalar form factor associated with the Faraday rotation term is given by 

eB 
m 



n p (k) 



1. 



where 



' B 



2.3 Axion Electrodynamics 



(2.6) 



(2.7) 



In this subsection we obtain an expression for the effective Lagrangian for axion-electrodynamics in a 
magnetized medium. It is important to note here that the axion contribution to the effective Lagrangian 
in a magnetized medium does get modified due to the presence of medium and magnetic field. However all 
the modifications can be included by redefining the axion-photon coupling constant [18]. The structure 
of the interaction Lagrangian remains identical to that at tree level. In light of this observation, we 
work with the tree level axion-photon Lagrangian in a magnetic field. In momentum space this effective 
Lagrangian is given by:, 



-Ank g^ 1 A v 



A^A V 



; T^k, t A v a 
1 2M a 



a(k 2 — m 2 )a 
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In the expression above, T^ v corresponds to the field strength tensor of the external field, is the vector 
potential for the photon, a is the axion field, 1/M a is the pseudoscalar-photon coupling and IP"' is the 
polarization tensor in matter along with the Faraday contribution, 

n^(fc) = U T (k)R^ + U L (k)Q^(k) + U p (k)P^ . (2.9) 

The equations of motion for the Lagrangian given by Eq. (|2.8p are the following: 

(-fc 2 .9 Q , + iWfc)) A"(k) = (2-10) 



(k 2 -m 2 )a^t b ^ A ^ k) (2.11) 
V ; 2M a V ; 

(2) 

From now on we would denote b a = k^T^a . 
2.4 Expanding A fl (k) in Orthogonal Basis. 

In this subsection we construct a system of orthonormal basis vectors out of the vectors available to us. 
Recall that the available vectors we have at our disposal are, it M the 4 velocity of the medium and fc M the 
external momentum of the photon. Using these two and T* 1 " we can further construct two more vectors, 

b {1 > = k^ v , (2.12) 

and 

bW v = k li P» u . (2.13) 

We also define, 

^-^rK- (2-14) 

We note that, the vectors b' 1 ', b^ and u are all orthogonal to k. Now we can construct another vector, 
I v , such that it is orthogonal to both u and b^ 2 \ 

F = (b&> - ^^«^ ■ (2.15) 

It is easy to verify that the vectors, l v u v , k v , w- x ' v are mutually orthogonal to each other. Therefore we 
can express the gauge potential as their linear combination, 

A a {k) =Ai(fc)Nib£> +A 2 (k)N 2 I a + A L (k)N L u a + k a N||A||(k). (2.16) 

The component, A\\(k) can be set to zero, since it is associated with the gauge degrees of freedom. In 
Eq. (|2.16|) the Nj's are normalization constants for the corresponding basis vectors. We can easily see 
that the normalization constants are given such to be, 

Nx = 1 (2.17) 

N2 = 7^ (2 - 18) 
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(2.19) 



The negative sign under the square root is for maintenance of the reality of the normalization constants. 
These constants take the values, 



Ni 
N 2 -- 



B Z K ± 

\K\ 

u)K±B z 



\K\ 



3 Equation of Motion in terms of the Form Factors 

The equation of motion for photon, Eq. (|2.10|) . can be written in a expanded form, 



[k 2 g^ - n T (k)R^ - U L {k)Q^{k) - U p (k)P„ v ] A u (k) = 



o \l a 
2M a 



(3.1) 



where we have used Eq. (|2.9[) . In order to arrive at the equations of motions in terms of the form factors, 
one needs to substitute Eq. (|2.16p in Eq. (|3.1|) and project out the different components of this equation. 
Multiplication from left by the normalized basis vectors, Ni x b\}\ N 2 x 1^ and Nl x leads to the 
equations, 



(fc 2 - U T {k))A 2 (k) + illpNtNi 



NiAi(fc) 



2M a 



(fc 2 - n r (fc))A!(fc) + iUpN^ 



(fc 2 -n L ) A L (k) 



A 2 (k)=0, 

iN L (bftu^ a 



2M„ 



respectively. The equation of motion for the pseudoscalar field can be expressed as, 



2M a 



N 2 A 2 (fc) 



2M a 



N L A L (k) 



2 2\ 

m a. 



(3.2) 



(3.3) 



This completes the closed set of equations involving the axion and the photon vector potential. It is 
easy to see that, if one sets the axion field to be equalto zero, one recovers the usual Maxwell Equations 
modified by the Faraday pieces. 



4 Mixing Matrix 

One can observe from Eqs. (|3.2|) . (|3.3|) that, the coupling to axions mixes the longitudinal component Al 
to the transverse components. So even if Al is zero to begin with, the same can be generated through 
coupling through the pseudoscalar field a. However the coupling is suppressed by the PQ symmetry 
breaking scale. One needs to study these equations carefully to determine the effect of pscudoscalar- 
photon coupling and magnetized medium on the Electro-Magnetic (EM) vector potentials. 

We are interested in a quasi-monochromative wave solution. We shall assume that the wave propagates 
in the z direction and express the solution in the form 

<fn(t,*) = e~ iut M0,*) (4- 1 ) 
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where (pi may represent the pseudoscalar field or any of the components of the electromagnetic wave. We 
work in the eikonal limit, where w « k is the largest energy scale. We may now express Eqs. (|3.2[) . (|3.3|) 
in real space in the matrix form 



[(- 2 + d 2 z )I M] 



( M{k) \ 
A 2 (k) 
A L (k) 
V a(k) ) 



= 0. 



(4.2) 



where I is a 4 x 4 identity matrix and the mixing matrix 

/ n T 



M 



\ 
















+n T 





4 2M„ 




.N 2 b^ ) I> J - 


n L 


.N L bfu" 
— 1 - 

2M Q 


2 M a 


2A/ a 



(4.3) 



/ 



In order to find out axion-photon oscillation, we need to diagonalize the matrix given in Eq. (|4.3p . 
Although it can be diagonalized exactly, the resulting formulae are too cumbersome to be directly useful. 
We make some simplifying assumptions to make the problem tractable. The longitudinal component has 
been shown to contribute negligibly in Ref. [51]. Hence at leading order we set Al(z) = and compute 
its contribution pcrturbatively. 

Our final aim is to compute the Stokes parameters I(z), Q{z), U(z), V(z) at some distance z, given 
the values of these parameters at the origin z — 0. For this purpose we may compute the density matrix 



' < A X A\ > < A 1 A* 2 > < A ia * > ' 
p(z) = ( < A 2 A\ > < A 2 A* > < A 2 a* > 
< aA* > < aA 2 > < aa* > 



(4.4) 



where the angular brackets <> represent ensemble averages. We work in the Lorentz gauge DA = and 
hence from these matrix we can directly compute the coherency matrix 



< ExEl > < E X E* 2 > 

< E 2 E\ > < E 2 E* > 



(4.5) 



where Ei are the electric fields. 



5 Solutions 

We now solve the equations of motion, ignoring the longitudinal mode. We work in the low temperature 
limit and ignore the terms proportional to temperature. We may express the resulting mixing matrix M 





f A 


iF 





M = j 


-iF 


A 


-iT 




K 


iT 


B 



where A = d*, B = m 2 , F = lulu cos 9 / (uj 2 — oj b ), and T = |i?|a;sin#/2M Q . Here 6 is the angle 
between the background magnetic field and the direction of propagation. Throughout we shall assume 
that the background magnetic field is independent of space and time. The component transverse to the 
direction of propagation is taken to point along the 'y' or '2' axis. 

We denote the eigenvalues of M by A,;. We define an alternate matrix M, 

I ix \ 

M = (M - AI)/(B - A) = I -ix -iy J , (5.2) 

V iy I ) 
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where I is an identity matrix, x — F/(B — A) and y = T / (B — A) . The eigenvalues Ai of M are related 
to Xi by 

A, = Xi(B - A) + A. (5.3) 
The eigenvectors of M (or M) may be expressed as 

/ x{Xi-l) \ 

|A»>=— -i\{Xi-l) , (5.4) 
1 \ vh I 

where D 2 = 2x 2 (Af — l) 2 + y 2 (2Xf — A,). The eigenvalues A; may be written as 



6> + 4tt\ 1 
3 




where 



Ai = 
A 2 = 

A 3 = 2v/-Qcos( ) (5.5) 



R 



i? = j-(9y 2 -18x 2 + 2) 
54 

^ 3 9 



(5.6) 




We point out that for a cubic equation we expect three real eigenvalues if Q 3 + R 2 < 0, which must of 
course be valid in our case. 

Using this we construct the unitary matrix U to diagonalize M and hence solve the propagation 
equation. The final result for the density matrix after propagating a distance z from the initial point is 

p{z) = UPU t p(0)UP' i U' ! , (5.7) 

where P is the propagation matrix, 

(exp(ikiz) 
exp{ik 2 z) | , (5.8) 



and ki = \Jlo 2 — Xi « u> — j 1 . 

Although the problem is solvable exactly, we find that there exist parameter ranges where the exact 
formulas are not very useful for numerical calculations. Furthermore in order to get some insight into the 
solution it is useful to explore it analytically in some limiting cases. The problem is a little tricky due to 
near degeneracy of the two modes of the electromagnetic wave. 

For astrophysical and cosmological applications we are usually interested in the case where both the 
parameters T and F are very small compared to \A — B\, i.e. x « 1 and y << 1. This is the case of 
weak magnetic field and pscudoscalar-photon coupling. In laboratory experiments the parameter y may 
not be small compared to unity if the plasma density is taken to be very small. We postpone a detailed 
discussion of laboratory experiments to a separate publication. For now we discuss the astrophysically 
interesting limit of x « 1 and y « 1. 
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We expand the eigenvalues in powers of x and y. This expansion is a little messy since we find that we 
need to include terms atleast up to order y 4 . This is because in the limit x — > 0, these terms contribute 

1 in the limit x 



at leading order. The expression Rj — Q 3 

R/ 

where 



1, y — > 1. Hence we may express 



-Q 3 = l 



0, 



(5.9) 



27x 2 /2 « 1 







27 



-.</ 



27 x 17 



54a; 2 y 2 << 1 



(5.10) 



The cos 1 (1 — a — 0) function has a non-analytic structure in the neighbourhood of unity. Expanding in 
powers of x and y we find [52] 

a + f3 3(a + /?) 2 



^(1 



0) 



V2 



G 



80 



(5.11) 



The eigenvalues are approximately given by 

2 



Ai 



^(x 2 + y 2 ) + ^= v / a + /3 + ^ 2 + y 2 f + 4f(^ 2 + + 



+ 



A 2 = 



A:< 



3V6 

^ + /3) + > 2 - 
2 V ^ 3^6 

1 / in a , 9 

27 (a + /?) + T8 ( ' T " 
1 + (x 2 + y 2 ) - l(a + 0)- \(x 2 + y 2 ) 2 - 1 -{a + (5){x 2 + y 2 ) . 



y 2 )\fu- 



(5.12) 



Here we have displayed terms accurate to order x 2 ,y 4 and x 2 y 2 . We again emphasize that in the limit 
x — ► 0, the dominant contribution comes from terms proportional to y A . Hence we need expand up to 
order y A . As we can see the expansion is quite cumbersome and does not yield a simple analytic result. 
However it is suitable for numerical calculations in the limit the parameters x and y are very small. We 
next consider some limiting cases in which a simple analytic result can be obtained. 



5.1 Limit 1: T « F « \A — B\ 

This is the regime where the Faraday effect is dominant and hence relevant at low frequencies. We expand 
the eigenvalues and eigenvectors in powers of T and obtain results accurate to order T 2 . This expansion 
may also be obtained directly by using perturbation theory. We may first exactly diagonalize the matrix 
M in the limit T = and then compute leading order corrections in T. The resulting eigenvalues are 
given by 

T 2 



Ai 



A — F 



A 2 = A + F 



2(A-F -B) 
T 2 



with the corresponding eigenvectors 



|Ai> = 



B 



1 

71 



2{A 
T 2 



F-B) 



T 2 



2{A-F-B) 2{A + F-B) 



(5.13) 



/ 1- 

i 1 - 

V 



4(A-F-B) 2 
A(A-F-B)' 2 



iF(A-F-B) 
' AF(A-F-B) 



A-F-B 
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|A 2 > = 



|As > 




2{A+F-B) 
T 2 

' 4(A+F-B) 2 



(5.14) 



4.(A-F-B) 2 

In the limit under consideration we may expand the denominators in powers of F/(A — B). Since 
F >> T, we compute the density matrix accurate to order T 2 /[F(A — B)] and drop terms proportional 
to T 2 /(A — B) 2 . We also assume that all correlators involving the pseudoscalar field are zero initially. 
The relevant elements of the density matrix are given by 



Pn(Z) 



o(Pn(0) 



P22(0)) 



cos(Ai 2 Z) 



Pia(O) 



--sin(A 12 Z) 



o(Pu(0) 



P12(0) 



sin(Ai 2 Z) 
2 

P12(0) 



P22(0)) 

sin(Ai 2 Z) 



4F(A - B 
cos(Ai 2 Z) 



(pii(0)-P2 2 (0)) 
iT 2 

■(1 - cos(A 12 Z)) 



+ c.c. 



(pn(0)-p 22 (0)) 



iT 2 



4F(A- 



(Pii(0)-paa(0)) 
- cos(Ai 2 Z) 



iT 2 



+ c.c. 



AF(A 
iT 2 



(1 - cos(Aia^)) 

(l-cos(A 12 Z))( jOll (0) + p 22 (0)) 



F(i4-B) 



sin(A 12 Z) 



Pi 2 (0) 



[-1 + cos(AiaZ)] (5.15) 



where Ai 2 = (Ai — A 2 )/(2w). We find the usual Faraday rotation along with additional contributions 
proportional to T 2 /[F(A — B)]. If the initial beam is unpolarized or linearly polarized, pseudoscalar- 
photon will generate circular polarization. The dominant contribution arises if the initial beam has a 
non-zero real part of pn or equivalently a non-zero Stokes parameter U . In this case in the limit Ai 2 Z << 
1, the circular polarization generated is proportional to T 2 Z/[u(A — B)]. Hence this contribution is 
proportional to w. 



5.2 Limit 2: F/\A-B\ «T 2 /\A-B\ 2 « 1 

We next consider another limiting case where the Faraday rotation effect is negligible and pseudoscalar- 
photon mixing dominates. This is the case where the leading order contribution leads to the standard 
pscudoscalar-photon mixing. In this limit we can also obtain the results by first diagonalizing the matrix 
M with F = and then obtaining leading order corrections in F . From the exact formulas we find that 
this expansion is valid only in the range, 



y 



F\A-B\ 

rp2 



« 1 



(5.16) 



It is clear that the range over which the leading order results are valid is considerably limited due to the 
nonlinear dependence on T in this inequality. 

The three eigenvalues, accurate to order (F/T) 2 are given by 



Ai 



A- 



A 



(A-B) 

T 2 
(A-B) 



F 
T 

-(A-B) 
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Aa = B — 



T 2 



(A-B) 



(5.17) 



We display the eigenvectors only to order Sf or upto terms linear in F, since these give the dominant 
corrections. The corresponding eigenvectors are given by 



|Ai > 



|A 2 > 



IA, > 



i8f 




1 

I>t~ 
1 

zv 



'S F 
-/' 
St 

St 



(5.18) 



where Dt = \Jl + S 2 ^, St = T/(A — B). The density matrix elements, accurate to order 8f are 



Pn{z) = pii(o) + (ipi 2 (o)<y F 

p 22 (0) 



1 -e 



i(Ai-A 3 )Z/(2u>) 



+ c.c 



P22{Z) 
P12(Z) 



D T 



1 + 4 e i(A 2 -A 3 )Z/(2 W ) 



+ (via(0)*ir [e^^"^)^^) - l] +c.c.) 



= -1(5f/3ii(0) 



1 _ e i(A 2 -Ai)Z/(2w) 



+ «<5f/O22(0) 



2 _ e i(A 2 -Ai)X/(2w) 



MO) 



e i(A 2 -Ai)Z/(2w) + £2 e i(X 3 -\ 1 )Z/(2u>) 



(5.19) 



5.3 Limit 3: F/\A-B\ ~T 2 /\A-B\ 2 « 1 



We finally consider the case where x and y 2 are of the same order but both are much smaller than unity. 
The eigenvalues A^ at leading order are 



Ai = ^x 2 + 2/74 - y 2 /2 
A 2 = -y/x 2 + yyA-y 2 /2 
A 3 = l + y 2 

with Ai related to A,; by Eq. 15.31 The corresponding eigenvectors, accurate to order y, are 

-1 



(5.20) 



[Ai > 



|A 2 > 



|A 3 > 



i 



— i\ 2 /x 
d2 1 



(5.21) 



where dx = yl + X 2 /x 2 , and e?2 = \/l + A 2 ,/^ 2 and c?3 = 1. The density matrix elements in this case are 
given by 



pu(Z) = pn(0) 1 



d\d 



[1 - cos(A 12 Z)] ) + p 22 (0)3232 [1 - cos(A 12 Z)] 



2 



11 



A 

xd\ 
2 



A 

xd\ 



xd\d\ 



cos(A 12 Z) - 



2y^ 



xd\d 2 



f^sin(A 12 Z) 



Pl2(0) + C.C. 



Pii( )^P| [1 - cos(A 12 Z)] + P22 (0) - _ [l _ cos(A 12 Z)] 

p i2 (0) + c.c 



• Al n 
• Xl n 



-iA 12 Z] 



4 



i 

4' 



,iA l2 Z 



1 

xdi 

A 

xd\ 



M 2 










\xd 2 j 


jLpiAj 


2 Z 


Pn(0) 



-iAi 2 Z 



xd2 

A 

xc? 2 



-iAi 2 Z 



3 -iAi 2 2 



p 22 (0) 

P12(0) 



[l-cos(A 12 Z)]pI a (0) 



(5.22) 



where Ai 2 = (Aj — A 2 )/(2cj). We point out that in the density matrix elements we have kept only 
the leading order terms. The terms of order y 2 have been dropped, unless they come multiplied by 
the distance factor Z . This is reasonable since for very large Z values such terms may give significant 
contribution but are otherwise expected to be negligible compared to the terms we have kept. 



6 Contribution of the Longitudinal Component 



In the full mixing matrix (|4.3p we have so far ignored the longitudinal mixing terms, i.e., M34 and M43. 
We now determine the contribution of the longitudinal mode by treating these terms as perturbation. 
The matrix element M34 of the mixing matrix (|4.3p is given by, 



|M; 



:S4 



2M 1 11 




The full 4x4 matrix, given in Eq. 14.31 may be expressed as, 



(6.1) 



M = 



M + M 
A iF 
-iF A 

iT 






n L 






-iT 









-iL 






iL 




(6.2) 



\B\ sin#/(2M) in the limit ui » uj v . Here we treat M' as a perturbation. The contribu- 



where L ks 

tion of the perturbation is small as long as 



<< 1 



(6.3) 



where A^ arc the unperturbed eigenvalues. In the limit of high frequencies, where the pscudoscalar-photon 
mixing effect may dominate, this difference is of order T 2 /A, ignoring the pseudoscalar mass. One can 
easily check that LA/T 2 << 1 for a wide range of parameters. In the low frequency limit, where Faraday 
effect dominates, the contribution of the longitudinal mode is also small. Hence we find that for a wide 
range of parameters we can treat the longitudinal mode perturbatively. 
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We next explicitly compute the modification of the eigenvectors and eigenvalues for the case discussed 
in section 5.3 above. Similar results apply in all cases. The eigenvectors and eigenvalues for the un- 
perturbed matrix Mq can be obtained from section 5.3. The unperturbed eigenvalues can be expressed 
as. 



A? = 



A 2 = 



2A 
T 2 

y2 

A 



A-\lF 2 
\lF 2 



2^4 
4A 2 



r 4 

AA 2 



(6.4) 



where we have set B = m 2 = 0. Similarly, it follows that, the unperturbed eigenvectors, 



IAS > 



|A^> 



1 



1 

~2 






\w 2 J 

Vi 





(6.5) 



where, v\ = iXi/x, v<i = iX^/x, 1)4 = —iy, w\ = yXi/x, W2 = yX 2 /x, W4 = 1, di = 1 and d\, d%, Ai and 
A2 are defined in section 5.3. We also have, 



A3 



|A^> = 





w 



(6.6) 



The leading order corrections to the eigenvalues are given by, 

M' jj = (x° j \M'\x° j ) 



(6.7) 



For all values for j, the leading order corrections to the eigenvalues due to the longitudinal pcrturbative 
piece vanish. The corrections to the eigenvectors are as follows, 



Cjk = 



A°|M'|A° 
(A fc -A 3 ) 



(6.8) 



We can verify that if k ^ 3 then the corrections to the eigenvectors are zero. Hence for j = 1,2,4, we 

have, 

|A, >= yfl - \C ]3 \ 2 |A« > +C j3 \X° 3 > . (6.9) 

where, 

„ —iLwj 
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Similarly, for | A3 > we can write, 

ia 3 >= A - E i c '3,i 2 iAg>+ J2 c ^ x °> ( 6 - n ) 

V J = {1,2,4} 3={1,2,4} 

We find that for a wide range of parameters, | Cj-3 1 = | C3 j- 1 << 1. Hence we see that the contribution of 
the longitudinal part is small and can be neglected at leading order. 

7 Conclusion 

In this paper we have given a general treatment of pseudoscalar-photon mixing in a magnetized medium. 
We solve the resulting coupled wave equations in several different regimes with the assumption that the 
frequency of the electromagnetic wave is much larger than the plasma frequency. The problem is a little 
complicated due to the presence of near degeneracy in the two transverse modes of the electromagnetic 
wave. We find, as expected, that at very high frequencies, the Faraday effect gives negligible contribution 
and the standard treatment of pseudoscalar-photon mixing is applicable. However at low frequencies, 
Faraday effect cannot be neglected. We extend the pseudoscalar-photon oscillation formulas so that they 
are applicable in this regime also. We find that in this case the oscillation effect is considerably modified 
due to the presence of near degeneracy in the mixing matrix. These results may be useful in many 
laboratory experiments and in astrophysical observations. 
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